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APPLICATION OF DIFFERENTIAL SIMILARITY TO FINDING NONDIMENSIONAL 
GROUPS IMPORTANT IN TESTS OF COOLED ENGINE COMPONENTS 

by James Sucec* 

Lewis Research Center 

SUMMARY 

The nondimensional groups required for similarity of flow and heat transfer in 
cooled engine components between a test rig and an engine are found by the application of 
the method of differential similarity to the governing partial differential equations and 
their boundary conditions for the hot gas, cooled metal, and coolant gas. The most im- 
portant of these groups required for similarity are geometric variables; Mach, Reynolds, 
and Prandtl numbers; the ratio of coolant to hot gas mass flow rate, the ratio of coolant 
to hot gas temperature, and the turbulence intensity of the free stream. Suitable re- 
arrangement of these groups makes it possible to arrive at the nondimensional groups 
normally used in tests of cooled engine components. The work in this report provides a 
sound theoretical basis for use of these groups, .which were selected in the past on the 
basis of heuristic, physical insight arguments and also because they are often the groups 
the experimenter has most control over. 

Differential similarity also uncovers some groups to which scant, if any, attention 
has been paid in the past. Apparently, these groups exercise a second order effect, at 
best. In addition, differential similarity is applied formally to two turbulent shear 
stress models, which causes the turbulence intensity of the free stream to surface as a 
nondimensional group. 


INTRODUCTION 

The work employs the technique of differential similarity to find the nondimensional 
groups which must have the same value for both a test rig and a cooled engine part in 
order that engine heat-transfer performance be properly modelled by a test rig. 

* Professor of Mechanical Engineering, University of Maine, Qrono, Maine; 
Summer Faculty Fellow at the Lewis Research Center in 1973. 



Often cooled engine components are tested at relatively low temperatures and pres- 
sures so that instrumentation and test equipment can be less complex and constructed 
less expensively than they would be for a test at the high-temperature, high-pressure 
engine conditions. Sometimes it is convenient to work with the full-scale hardware, for 
instance, when one is interested in the temperatures in a first stage vane. At other 
times a larger than prototype size scale model may be dictated by a desire to investigate 
more detailed, basic phenomena such as the influence of a pressure gradient on the film 
effectiveness very near a film slot on a surface. In both cases one must be able to de- 
sign the test rig properly and then scale the experimental information collected from the 
test so as to predict prototype behavior. This requires that the test rig (model) be sim- 
ilar to the engine component (prototype). 

One necessary condition for this similarity is geometric similarity, that is, the use 
of a scale model. The other conditions for similarity can be found in a number of ways 
including the use of the Buckingham pi theorem, examination of available exact analytical 
solutions to closely related problems, use of physical insight and reasoning, and exam- 
ination of the governing partial differential equations and the boundary conditions. This 
last technique can be employed even though the equations and boundary conditions may 
not be presently amenable to an exact analytical solution. In reference 1 a combination 
of physical reasoning and some nondimensional groups presented in reference 2 leads to 
some of the conditions necessary for similarity. The groups used and the scaling pro- 
cedure recommended are checked by comparing measured airfoil temperatures at high 
gas temperature and pressure with scaled up values based on tests at low gas temper- 
ature and pressure. The agreement seems reasonable in the range of conditions ex- 
plored. Chapter 11 of reference 3 presents similarity conditions based on heuristic ar- 
guments, physical insight, and the author's experience and knowledge of what is impor- 
tant when investigating forced convection heat-transfer phenomena. Two tables of sim- 
ilarity states valid for a high-performance jet engine at takeoff conditions and at cruise 
conditions are presented. Reference 4 is a much more detailed and expanded version of 
the similarity material presented in reference 3. The basis for selection of the appro- 
priate nondimensional groups is still the same as in reference 3, and the steps one would 
follow in using them are thoroughly explained. In addition, the authors emphasize some 
of the pitfalls that await those who violate the major constraints imposed by similarity 
requirements. The authors point out the serious errors that can occur if one fallaciously 
reasons that cooling performance can be demonstrated by a test at near prototype tem- 
peratures but at low pressures very near the standard atmosphere. The procedures and 
the caveats advanced in reference 4 constitute an invaluable addition to the literature for 
those involved in testing the cooled components of the engine. However, there is still a 
gap in the general area of a rigorous foundation for the nondimensional groups which are 
sufficient for similarity. 
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This report uses a rigorous mathematical method, differential similarity, on an 
idealized but still fairly general physical model to find the nondimensional groups re- 
quired for similarity. Once these groups are found, they are examined to see if they 
imply the groups found by the less rigorous approach of reference 4; since they do, they 
can serve as the underlying framework of the structure presented there. Additionally, 
differential similarity is formally and explicitly applied to some turbulent shear stress 
models to see if any additional groups appear over and above those which exist in a 
purely laminar flow. 


ANALYSIS OF STEADY, LAMINAR FLOW 

To find the nondimensional groups which ensure similarity by the use of differential 
similarity, one writes down all the partial differential equations, the algebraic equations, 
and the side conditions that govern the solution of the problem at hand; makes all equa- 
tions and side conditions nondimensional; and then by direct inspection writes down the 
implicit functional dependence of all the dependent nondimensional groups upon the inde- 
pendent groups and the other nondimensional parameters. For- more 'detail, one may 
wish to consult references 5 and 6. 

Although the results of the following analysis are valid for most air-cooled station- 
ary components of an engine, the following idealized model is considered for concrete- 
ness. A hot gas flows in a steady, laminar fashion through a cascade of stationary metal 
vanes inside of which a coolant gas flows in a steady and laminar mode. Body forces are 
neglected, and both the gas and coolant are assumed to obey the perfect gas equation of 
state, have internal energy per unit mass and enthalpy per unit mass which depend only 
on temperature, have a dynamic viscosity which is only temperature dependent, and be- 
have as linearly viscous fluids in regard to their viscous stresses. For simplicity, the 
thermal conductivity, dynamic viscosity, and specific heat of both the gas and the coolant 
have their temperature dependence given by power law forms. The metal is considered 
rigid and isotropic, and a power law form is used for the temperature dependence of its 
thermal conductivity. Gas radiation effects and mass diffusional effects are not dealt 
with. These constraints define the idealized model which is now rigorously analyzed for 
the nondimensional groups which ensure similarity between engine conditions and test 
conditions. 

Later the generalization to steady turbulent flow is made, and this requires addi- 
tional assumptions to be invoked because of the need for a choice of turbulence model. 

A schematic of one of the vanes under consideration and some of the associated 
nomenclature are presented in figure 1. (All symbols are defined in appendix A. ) 

The governing equations for the flowing hot gas are as follows: 
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Conservation of mass: 


9(pu) , 3(pv) H 3(pw) 
ax ay az 


(i) 


x-direction Navier -Stokes (see ref. 7): 
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(The Navier -Stokes equations for the y- and z-directions are similar in form to eq. (2) 
and are not shown since they do not yield any nondimensional groups that are different 
from those which can be extracted from eq. (2). ) 

Thermal energy: 



ax ay az 

(3) 
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State: 
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(4) 

(5) 
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It is assumed that the temperature dependent fluid properties, dynamic viscosity p, 
thermal conductivity k, and specific heat at constant pressure c p can be adequately 
represented in the temperature range of interest by the following power law forms: 


M(T) = fJ-Q 



( 6 ) 


k(T) = k Q 



(V 


c p (T) = c 


p,oo 



For the metal, which is assumed isotropic and rigid, one has the following 
equations: 



( 8 ) 


( 9 ) 



( 10 ) 


For the coolant one has equations that are the same as equations (1) to (5) with all quan- 
tities except the independent variables x, y, and z having a subscript c to identify 
them as being local values within the coolant. Hence, the equations for the coolant are 
not shown explicitly. Also, for the coolant we have 



(ID 



( 12 ) 
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c p,c c p,°°,c 


(13) 



At this point it might be noted that the coordinate system shown in figure 1 is rec- 
tangular Cartesian, so that the governing partial differential equations are relatively 
simple and have no curvature terms. In the event that the circum stances allow thin- 
boundary-layer flow, the coordinate system can be viewed as a curvilinear one which 
follows the vane surface. 

The boundary conditions must be written down. For the hot gas as x — -°° 


U =U oo P=P oo' 

v = 0 T = > 

w = 0 


(14) 


On all solid surfaces 


u = v = w = 0 


(15) 


and downstream in the gas as x 


3u _ 3v _ 3w _ o 
9x 3y 3z 


3T 

3x 


= 0 


(16) 


(Note that a precise specification of the downstream boundary, x — is not important 
unless the velocity and temperature fields are strongly elliptic in the x-direction. If the 
flow is of the thin -boundary -layer type, nothing need be specified as x — °°, since the 
equations are then parabolic in x. Hence, in eq. (16) relatively weak conditions are 
used on the downstream boundary. ) 

For the coolant the boundary conditions are 
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z = 0 


On all solid surfaces 


u = v 
c v c 


= 0 


w c = w ci 



(17) 


u 


c 



(18) 


The remarks previously made in connection with the downstream boundary for the hot 
gas are taken as applicable to the coolant as well. 

At the interface between the gas and the solid metal the conjugation conditions are 


T 



(19) 


- k m — = " k — + e m a 
m Sn 9n 1 



( 20 ) 


The second term on the right side of equation (20) represents the net radiant loss per 
unit area by the metal due to gray body radiant exchange with solid surfaces. In this 
term B„ . is Gebhart's absorption factor, the fraction of the radiant energy emitted by 
the surface m which is ultimately absorbed by the j surface. 

Two conjugation conditions of the form of equations (19) and (20) also hold for the in- 
terface between the coolant gas and the metal. 

Equations (1) to (20) are made nondimensional by dividing all velocities, whether in 
hot gas or coolant, by u M (hot gas velocity upstream of cooled component); dividing all 
temperatures, hot gas, coolant gas, and metal, by Tq (inlet stagnation temperature of 
hot gas); and dividing all space coordinates by the characteristic length L, which is 
chosen as the airfoil chord in this report. Also, pressure and density are made dimen- 
sionless through division by P M and p m , respectively, the static conditions at the hot 
gas inlet. Differential similarity methods next require formal transformation of the di- 
mensional dependent and independent variables to the nondimensional ones in all the equa- 
tions. Once this is complete, the equations are made nondimensional. A representative 
illustration of this procedure is presented in appendix B. 

When equations (6) to (8) are combined with equations (2) and (3), equations (11) to 
(13) are combined with the counterparts of equations (2) and (3) for the coolant, and 


7 



equation (10) is inserted into (9), the method of differential similarity finally yields the 
following nondimensional equations: 

For the hot gas : 

Conservation of mass: 


3(pU) + 9 (pV) + 3(pW) = 0 

ax aY az 


( 21 ) 


x -direction Navier -Stokes: 


P 






0 d / 3U + other velocity \ 

/ ax Re ax 

V 0X gradient terms/ 


( 22 ) 


Thermal energy: 



+ EuEc 


U^ + 

\ ax 


Re 



(23) 


For the sake of compactness of form, terms with the same general form and char- 
acteristics, such as the three terms in parentheses on the left side of equation (22), are 
often represented by + . . . This, of course, is no loss once all the nondimensional de- 
pendent and independent variables have appeared, since these additional terms do not 
contribute any new nondimensional groups. 

State: 


p = - EuEc 
0 



(24) 


For the metal: 


_3_ 

ax 



= o 


(25) 
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For the coolant: 
Conservation of mass: 


±(p c U )+... =0 

ax c c 


(26) 


x-direction Navier -Stokes: 


- / 0U c 


3P C M c o 1 1 

-Eu - + c > u _L _L 

0X ^Xq Rb 0X 


d c / 3U 

e c 2 —1 + 

\ ax 


(27) 


Thermal energy: 
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(28) 
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C 0 c R c T 0 


(29) 


The nondimensional boundary conditions for the hot gas are 
As X - -oo 


U = 1 


V= 0 


W = 0 


As X-oo 


On solid surfaces 


P = 1 


e = 


au _ dv aw = _ae = 0 
ax ay ~ ez ax 


u = V = W = 0 


(30) 
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In the coolant the boundary conditions are 
At Z = 0 


On solid surfaces 


w 


W c = 


Cl 


U 



v c =o 
u c = 0 



> 



( 31 ) 


Conjugation conditions on solid surfaces are 
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Inspection of equations (21) to (32) reveals that the nondimensional dependent varia- 
bles are a function of the nondimensional independent variables and all nondimensional 
parameters. Hence, 
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(u, v, w, P, e, P , e m ,v c , v c , w c , P c , e c , p c ) = G, 
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( 33 ) 


where 


Eu = 


P 00^00 


Pr = 


_ ^0 c p, °° 

k 0 


Re = 


Poo U cc L 


Mq 


U 


Ec = 


C p, «> T 0 


The function G i in equation (33) is different for each of the dependent nondimen- 
sional variables on the left side of the equation. 

Basically equation (33) means that the nondime nsional fields of velocity, temper- 
ature, pressure, andj density in the hot gas and in the coolant, as well as the nondimen- 
sional metal temperature field, are the same in both model and prototype, test rig and 
engine, if the corresponding nondimensional groups in the test rig and the engine which 
appear in the argument of the G i function are the same in the test rig and the engine. 

In the framework of the problem addressed in this report, this requirement and geo- 
metric similarity constitute sufficient conditions for similarity between the test rig and 
the engine. 


EXTENSION TO STEADY ON THE AVERAGE TURBULENT FLOW 

As is well known, for a steady turbulent flow, equation (2), for instance, contains 
additional momentum flux terms due to the fluctuation components of velocity. It is 
usual to treat these terms as additional apparent stresses, called turbulent, or Reynolds, 
stresses. To see whether differential similarity will predict any additional nondimen- 
sional groups due to turbulence, consider, for the sake of simplicity, the case of 
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incompressible, thin -boundary-layer flow. Calling the apparent turbulent shearing 
stress in the x-direction x , one can show that equation (2) reduces to 


P 






ay 



( 34 ) 


As shown in appendix B, for this equation without the turbulent stresses, the non- 
dimensional form becomes 


U*U. 

ax 


= -Eu 


d -* + ±±(> 

3X Re 3Y \ 


aY 


a t. 


2 

oo oo 


p,- u ' 


t,x 

sy 


(35) 


Now if one chooses to use mixing length theory together with an algebraic prescrip- 
tion of the mixing lengths for a boundary layer near a single wall, one obtains the follow- 
ing equations, according to reference 8, page 29: 


Since 


T tx =p„(0.41)V 

/ 3u\ 2 
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0 < y < °* 09 6 
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^-pJD.OOjV 
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U=u/u oo and Y = y/L and A can be defined as 6/L, equations (36) become 


V .«4<M«V(g)[ 


T t,x = ^ 4 < 0 -° 9 > 2 A 2 (^) 2 


2 


(37) 


Differentiating equations (37) with respect to y and then changing to the nondimensional 
coordinate Y give 
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( 38 ) 


!^ = ^<o.4i) 2 ±[y 2 (wf| 
ay L 3Y L VaY / J 

► 

!V = ^(0.09) 2 A 2 -l[('wfl 

ay l 3Y [VaY / J J 

2 

When equation (38) is inserted into (35), the p^u^/L term cancels, and the nondi- 
mensional equation, including the turbulent stresses, has one of the following two forms, 
depending on the position within the boundary layer : 


u™ + .. 
ax 


= -Eu + 

ax 


±±(e d ™ + ...\ 
Re 3Y V SY J 



(0.09) 2 A 2 _L / + . .. 
[[ dY\dYJ 




> 


J 


(39) 


Upon examination of equation (39), one sees that no new groups are introduced by 
taking explicit account of the turbulent stresses, at least not for the simple turbulence 
model involving prescription of the mixing length through algebraic equations. (Since A 
is defined as the nondimensional thickness of the boundary layer, it depends on some of 
the same groups as U and hence introduces no new groups. ) 

With regard to models for the turbulent shearing stress, the next level of sophisti- 
cation is the so-called one equation model of turbulence, where a length scale of turbu- 
lence and a characteristic scale of the turbulent velocities are used. The length scale is 
taken to be the usual algebraic mixing length prescription, while the velocity scale used 

is proportional to the square root of the turbulent kinetic energy. When K is defined as 

2 

the average turbulent kinetic energy divided by u w , the nondimensional form of the par- 
tial differential equation describing the distribution of K can easily be shown to be 


«« + ... = ± fey* fi as) + Lmlx ft (ssf - CpK3/2 

3X 8Y\ ® k 8Y j m “T \8Y/ L mix 


(40) 


(See ref. 8 for the dimensional version of the turbulent kinetic energy equation. ) 
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The boundary conditions in nondimensional form for this equation are 
On all solid surfaces "'i 


K = 0 


As X — °° 


As X - -<*> 


9K = 0 

0X 


( 41 ) 


K 


= K 


j 


In equation (40), L mix is the nondimensional mixing length, and reference to equa- 
tion (39) and the immediately preceeding work indicates that the presence of L mix does 
not give rise to any new groups. However, equation (40) gives one new dependent group, 
K, to be added to the groups on the left side of equation (33). The C D in equation (40) 
is a pure number, while cr k is a sort of turbulent Prandtl number for turbulent kinetic 
energy. It is the ratio of the eddy diffusivity of momentum to the eddy diffusivity of tur- 
bulent kinetic energy. This cr k and the ordinary turbulent Prandtl number a t are two 
additional groups that should be added to the right side of equation (33); <x t appears in 
the nondimensional thermal energy equation for turbulent flow when one relates the tur - 
bulent transport coefficients of heat and momentum. Also, the nondimensional boundary 
conditions lead to a new group to be added to the right side, namely, K m . But is a 
multiple of the more common turbulence parameter, the turbulence intensity TI^, so 
this is used in its place. 

In appendix C it is shown that Tq/T^, Eu, and Ec are functions of inlet gas Mach 

number M , of l, and of c /R, which itself depends on the specific heat ratio of the 

P J 

hot inlet gas. Utilizing these results and the ones just derived for turbulent flow, one 
obtains 


14 



(u, v, w, p, 0, , U c , V c , w c , P c , e c , 2 C , K) = 


y y 7 . Pi* m ^ i <r C P>°° m0 ^c,0 ^ ? c,0 _ c p,°°_ 

A ) i )^ J 5 IVC ) ri ) lvl ooJ U J £ ')b) _ J u > “ 5 J 

0 c p,«,c 


n> 

H kp iiq k n c. 


„ R w ci T cl P ei £ m° rT 0 £ m" T (L k mO £ m <7T 0 

c ’ R? 77’ 77 P? ~kT’ ~T~ 7“’ T7’ 


c U co A 0 00 0 

fe m aTn \ 

"k^oo 

c ,o r 


v 0 


c,0 c, 0 


(42) 


The function Fj is different for each of the nondimensional dependent variables. 


RESULTS AND DISCUSSION 

One of the major results of this analysis is equation (42); for a test rig and an en- 
gine that are geometrically similar this equation gives the sufficient conditions for equal- 
ity of the nondimensional temperature, pressure, density, and velocity fields in the hot 
gas, the metal, and the coolant. These sufficient conditions for similarity between the 
test rig and the engine are the equality, in the test rig and in the engine, of every nondi- 
mensional group in the argument of the functions on the right side of equation (42). As 
is generally true, one can have strict similarity, that is, meet all the conditions for 
similarity, only if the test rig is the actual engine operating under actual engine condi- 
tions. However, as mentioned previously, many of these requirements exert only a 
slight influence on the dependent variables and sometimes may safely be dispensed with. 
Brief arguments are now presented for the dropping of many of these constraints. In 
this discussion the superscript t refers to the test rig (model), and the superscript e 
refers to the engine (prototype). 

Perhaps the most convincing argument that many of the groups on the right side of 
equation (42) exert only a weak influence on similarity between the test rig and the en- 
gine, and on proper scaling, is that tests have been run without equality between the test 
rig and the engine for many of the groups and that direct experiment has shown that the 
lower pressure and temperature data do scale up and agree well with the measured 
values at engine conditions. This is shown in reference 1. 

One of the requirements for strict similarity from equation (42) is that 

d* = d e (43) 
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That is, the exponent on the power law fit to the temperature dependent dynamic vis- 
cosity of the hot gas must be the same in the test rig as in the engine. A plot of the dy- 
namic viscosity of air from the values given in reference 5 yields the conclusion that a 
single choice of d in equation (6) adequately represents p. over a wide temperature 
range. Hence, use of air in both the test rig and the engine tends to guarantee at least 
approximate equality of the two members in equation (43). If air is also the coolant, it 
seems reasonable to suppose that the nondimensional groups d, l , g, Cp ^/R, ju c q/hq, 
d, Z„, k n An> «/ c ™ 311(1 R/R„ would have the same value in both the test 

rig and the engine. 

Another ameliorative influence is the fact that satisfying the most important condi- 
tions for similarity ensures the same temperature ratio across the boundary layer at 
corresponding points in the test rig and the engine. This condition gives the correct ef- 
fect of temperature dependent properties on both the surface coefficient of heat transfer 
and the wall shearing stress for a large number of different gases, as is pointed out in 
reference 7. The available experimental evidence and the present thinking (ref. 8) indi- 
cate that cr^ and are practically equal and do not vary much from fluid to fluid ex- 
cept for the extreme Prandtl number fluids such as liquid metals and hydrocarbon oils, 
neither of which would be used as the test fluid for the engine hot gas or the engine cool- 
ant. Hence, the equalities that follow are automatically satisfied: 



It is usual to select test fluids and materials that are close enough to engine condi- 
tions to allow the following equalities: 



(45a) 


b 


t 


= b 


(45b) 



(45c) 


Actually equation (45c) is automatically satisfied as long as equation (45a) is satisfied in 
conjunction with a relation previously discussed, namely, 
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Next consider the groups which would ensure similarity of the solid to solid radia- 
tion, the groups in equation (42) which contain the total hemispherical emissivity of the 
metal e ffi . The absorption factor B m j depends on angle factors and the total hemi- 
spherical emissivities of all the solid surfaces which interact radiantly. Because of 
geometric similarity, corresponding angle factors in the test rig and in the engine are 
equal. If the same solid materials are used in both the test rig and the engine, approxi- 
mate equality of corresponding emissivities is obtained unless the emissivities are very 
temperature dependent. In that case, the lower solid surface temperatures of the test 
rig as compared with the engine make equality of the emissivities difficult to obtain. 
However, these problems pale in comparison with the difficulties raised by the presence 
of Tq in each of the radiation terms. For instance, strict similarity requires that 



Cancelling out the Stefan -Boltzmann constant a and rearranging give 


(45d) 



(45e) 


Despite the moderating influence of the thermal conductivity ratio, the temperature 
ratio to the fourth power dominates the right side of equation (45e) whenever the hot gas 
test rig temperature is markedly lower than the hot gas engine temperature. Obviously, 
equation (45e) can then easily require values of greater than 1 in order that equa- 
tion (45d) be satisfied, but such values are thermodynamically forbidden. The diffi- 
culties involved in satisfying the radiation requirements for similarity are pointed out in 
reference 4. 

Dropping all the radiation groups, as well as the other groups previously discussed, 
simplifies equation (42) and yields 
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(46) 


(U, v, w, p, 6 , 5 , e m , U c , V c , W c , P c , e c , p c , K) = 


w . T . P . 

= Pl(X,Y,Z,Re,Pr,M M , -Hi, _2, TI^ 
' u » T o P » 


Thus, when air or gases like air are used for the test fluids and the same metal as in the 
engine is used in the test rig, approximate similarity between the test rig and the engine 
is enforced by equality of the nondimen sional groups on the right side of equation (46) for 
the test rig and the engine. 


Illustration of Scaling Procedure for Hot Gas Side Surface Coefficient of Heat Transfer 


To understand how the method of differential similarity is used to predict engine per- 
formance based on available experimental results from the test rig, consider the prob- 
lem of finding the engine gas side surface coefficient h e when the measured values from 
the test h* are available. One has, from the definition of h, 


h = 



(47) 


where T m is the local metal temperature, and T g is a static temperature outside the 
boundary layer, or a local bulk mean temperature if the flow is not of the boundary -layer 
type, or an aerodynamic recovery temperature. This last choice is the one ordinarily 
made since it allows the high-speed h to be practically the same as the low-speed h. 
The precise choice of T_ is relatively unimportant to the discussion that follows as long 
as one is consistent when using the result. Utilizing equation (47) along with the defini- 
tions of N q and e causes (47) to become 


h 


L l 0N„ 


N o=° 


e m ” 6 s 


(48) 


or after rearrangement 
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( 49 ) 




e 

m 



K the groups on the right side of equation (46) are equal in the test rig and the en- 
gine, it then follows from equation (46) that e m , 6 S , and (30/3 N q ) n _q are exactly the 

same functions of the nondimensional space coordinates X, Y, and Z in both the test 
rig and the engine. Therefore, hL/k^ has the same value at corresponding points for 
the test rig and the engine, and the following form is obtained for scaling purposes: 


h e L e 


hV 


(50) 


Equation (50) is used to predict h e once h* is measured. A similar procedure can be 
employed to scale coolant side conditions such as the surface coefficient and the wall 
shearing stress. Because of many analytical and experimental results (see refs. 6 and 
7), it is well known that the local Nusselt number depends upon the molecular Prandtl 
number to about the one -third power, at least in the moderate Prandtl number range. 
This information can be used to write down the form equation (50) takes if strict Prandtl 
number similarity is not achieved: 


^ = 5V:/'2£f\ 1/3 (si) 

k ? k o W/ 

Equation (51) is essentially equivalent to equation (11) of reference 4, which evolves from 
physical insight arguments. 


Reduction of Groups to More Useful Forms and Comparison 
with Groups from References 

Some of the groups on the right side of equation (46) are not in the form which is 
most easily controlled by the experimenter when he sets the test conditions that pre- 
sumably effect similarity between the test rig and the engine. The tasks of reduction of 
groups to more useful forms, comparison with groups used by other investigators, and 
verification of the groups arrived at by other investigators are addressed next. 
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From equation (46) one of the requirements for similarity is that 


or 


Another is that 


These lead to 


Re e = Re 1 


u e L e 

OO OO 


e 

^0 


t t t t 

Poo U oo L 


t 

M 0 





u: 


c 


t 

oo 


The perfect gas equation of state and equation (53) inserted into (52) yield 

R e p L <£ L ‘ t' mo 

Using the relation for c x and T^/T^ developed in appendix C, the fact that 



from equation (42), and the result that 



(which can be shown) enables one to rewrite equation (54) as 


(52) 


(53) 


(54) 


(55) 


(56) 
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(57) 


_0 Mq 

3 e t 
0 M 0 



However, since a power law dependence of viscosity upon temperature was postulated, it 
can be easily shown that 


e 

Mo 

t 

MO 


e 

Moo 

~7 

Moo 


(58) 


Combining equations (57) and (58) yields 

(59) 

PS t i V **0 

Equation (59) and its equivalent (57) are consequences of the equalities (52), (53), and 
(55) predicted by the method of differential similarity. Equation (57), in particular, is 
useful for finding the required Pq for similarity. Ordinarily one would know Pq, Tq, 
R e , R*, and lVl 6 ; would choose a Tq less than T^; would determine /j.q and ji^ 
from property tables; and would then solve equation (57) for Pq, a quantity easily con- 
trolled and set by the investigator in the test rig. 

Comparing equation (59) with equation (5) of reference 4, one sees that they are iden- 
tical once lVl 6 is set equal to unity (because the ref. 4 equation is based on full-scale 
hardware) and once rfyr e in reference 4 is set equal to unity. The latter ratio is a 
function of the specific heat ratios in test rigs and engines, and the similarity require- 
ment of equation (55) forces the specific heat ratio to be unity. Equation (5) of refer- 
ence 4 is derived by the approach previously referred to as a heuristic, physical insight 
combined with experience and knowledge type of approach. The differential similarity 
analysis which led to equation (59) serves as a more rigorous basis for and substantia- 
tion of the physical approach. 

Using expressions for the various mass flow rates, the perfect gas equation of state, 
the similarity requirements of equation (46) or (42), particularly the requirements 





u 


(60) 
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( 61 ) 


and 


one can show that there results 





(62) 


(63) 


Thus, the ratio of coolant mass flow rate to hot gas mass flow rate in the test rig 
must equal the mass flow rate ratio in the engine. This consequence of the differential 
similarity analysis is often used as a starting point in some of the more physical ap- 
proaches to similarity as, for instance, in reference 1. Equation (63) also agrees ex- 
actly with its equivalent in reference 4 , which is derived there as a byproduct of another 
physical requirement. 

One can also show from the similarity requirements, equation (46) or (42), that at 
corresponding points in the flow (such as a film cooling slot, for instance) the momentum 
thickness Reynolds numbers, momentum flux ratios, and mass flux ratios are the same 
in the test rig as in the engine. That is, 

Re 6. = Re s. (64) 



(65) 


( 66 ) 


where |v| is the magnitude of the nondimensional velocity vector V, and terms like 
are to be interpreted as the vector dot product, = V • V, not as the square of the non- 
dimensional Y component of this vector. Requirements (64), (65), and (66) along with 
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geometric similarity assure the same adiabatic film effectiveness in the test rig and in 
the engine. By way of contrast, reference 4 asserts it is self-evident that equalities 
(64), (65), and (66) are physically reasonable requirements to force similarity. Finally, 
it can be shown that the similarity requirements (eq. (46) or (42)) lead to equation (9) of 
reference 4 and to the result that cp (a commonly used dimensionless wall temperature) 
is the same in both the test rig and the engine: 

< p e = (67) 

Hence, the heuristic, physical insight approach to similarity constraints as ad- 
vanced in reference 4 is totally verified by the present work, which represents a more 
rigorous, more mathematical approach to the similarity problem. The present work 
predicts the need to maintain the same turbulence intensity in both the test rig and the 
engine, as does reference 4 also. In addition, the method of differential similarity pre- 
dicts numerous other groups, on the right side of equation (42), that theoretically should 
be held constant between the test rig and the engine. These other groups, some of which 
have not appeared in the other references, most likely are of second order effect in most 
instances for the reasons noted previously. They could, however, assume a greater im- 
portance if any drastic changes were made in the types of fluids and metals used in the 
test rig relative to the types used in the engine. 

Equations (57) and (62) are two of the most useful results of this report for the de- 
termination of what pressure and temperature levels to set in the test rig to achieve sim- 
ilarity with the engine. However, since these equations are equivalent to the relations in 
reference 4, it is recommended that the reader follow the procedure given in reference 4 
to set the test conditions. Representative numerical calculations employing equations 
(57) and (62) have been made by the author, and agreement is good when compared with 
the similarity states given in table II of reference 4. The very slight differences can be 
traced to the slightly different viscosity -temperature relations used and the fact that ref- 
erence 4 makes a correction for the small difference in specific heat ratios between test 
rig and engine conditions whereas the method presented in this report does not. How- 
ever, the correction for different specific heat ratios could be incorporated easily, since 
it did originally appear in the derivations, but was cancelled out by the similarity re- 
quirement (55). 


SUMMARY OF RESULTS 

The method of differential similarity was used to rigorously derive nondimensional 
groups that are sufficient to ensure that the nondimensional fields of temperature, ve- 
locity, pressure, and density are the same in both a test rig and an engine for the hot 


23 



gas, coolant, and metal components. The following results were obtained: 

1. Arguments were given based on available experimental evidence and on the usual 
choice of types of fluids and metal used; these arguments reduce the large number of 
nondimensional parameters sufficient for similarity to a much smaller necessary num- 
ber of groups for approximate similarity in many cases. 

2. It was shown that approximate similarity, both aerodynamic and thermal, can be 
achieved in a similar geometry if the inlet Mach number, Reynolds number, ratio of 
coolant temperature to hot gas stagnation temperature, ratio of coolant to hot gas mass 
flow rate, Prandtl number, and free -stream turbulence intensity in the test rig are set 
equal to their counterparts in the engine. 

3. The derivation and the results presented serve as a proof of the validity of the 

similarity groups proposed by earlier investigators using heuristic, physical insight 
arguments. .... 

4. The method when applied to some simple turbulence models predicts additional 
groups needed for similarity in steady turbulent flow, although only one of them, the tur- 
bulence intensity, is thought to be important to the usual type of experiment. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 2, 1976, 

505-04. 
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APPENDIX A 


SYMBOLS 


m,j 


b 

C D 

c 

C P 

d 

Ec 

Eu 

F j 


f l’ f 2’ f 3 


G i 

g 

h 

i 

i 0 

j 

K 

k 

L 

L 

l 

M 

m 


mix 


n 


N„ 


absorption factor 
exponent 

constant related to dissipation of turbulent kinetic energy 

local speed of sound 

specific heat at constant pressure 

exponent 

Eckert number 

Euler number 

function 

functions 

function 

exponent 

surface coefficient of heat transfer 
static enthalpy per unit mass 
stagnation enthalpy per unit mass 
index 

nondimensional time-averaged turbulent kinetic energy 

thermal conductivity 

airfoil chord 

mixing length 

exponent 

Mach number 

nondimensional space coordinate measured perpendicularly outward from 
metal surface into coolant 

space coordinate measured perpendicularly outward from metal surface 
into hot gas 

nondimensional n 



P local pressure 

Pj function 

a 

P nondimensional pressure 

Pr Prandtl number 

p number of radiation surfaces on hot gas side 

P c number of radiation surfaces on coolant side 

R particular gas constant 

Re Reynolds number 

5 entropy per unit mass 

T static temperature (absolute) 

Tq stagnation temperature (absolute) of hot gas at inlet 
TI^ turbulence intensity at inlet 

U nondimensional x-component of velocity 

u x-component of velocity 

V nondimensional y -component of velocity 

v y-component of velocity 

W nondimensional z -component of velocity 

w z-component of velocity 

X nondimensional x-coordinate 

x space coordinate 

Y nondimensional y-coordinate 

y space coordinate 

Z nondimensional z -coordinate 

z space coordinate 

T function of y 

y specific heat ratio 

A nondimensional velocity boundary -layer thickness 

6 local thickness of velocity boundary layer 
local momentum thickness 

e total hemispherical emissivity 
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e 


temperature ratio, T/Tq 
p. dynamic viscosity 

£ dummy variable for T 

p mass density 

ft nondimensional density 

cr Stefan -Boltzmann constant 

cr ^ effective Prandti number for diffusion of turbulent kinetic energy 

CT t turbulent Prandti number 

r shearing stress 

$ dissipation function 

w mass flow rate 

Subscripts: 

c coolant 

ci coolant inlet 

i index 

m metal 

mO metal condition at hot gas inlet stagnation temperature 
ref reference state 

s undisturbed conditions far from surface 

t turbulent 

x x -direction 

0 stagnation conditions at hot gas inlet 

<» static condition of hot gas at inlet 
Superscripts: 
e engine (prototype) 

test or test rig (model) 


t 



APPENDIX B 


ILLUSTRATION OF DIFFERENTIAL SIMILARITY 
APPLIED TO MOMENTUM EQUATION 

Equation (2) is the x-direction Navier -Stokes equation for the hot gas, and equa- 
tion (6) is the viscosity -temperature relation for the hot gas. When these equations are 
combined, and only one of the viscous stresses is retained explicitly (because the others 
do not add any new nondimensional groups), the result is 


/ 3u 

011 

3u N 

i s _ 0P + _0_ 


/ T 

\ d 9u + 

u — + v — + w 



Mq 


i — + . . . 

\ 3x 

3y 

3z / 

3x 0y 

\ T o 

/ ay 


Equation (Bl) is transformed to nondimensional dependent variables by introducing 
the following definitions: 


e =-l 

T o 


W=?L 

U . 


U = 


u 

u„ 


p = 


V = — 
u . 


~ p 
p = — 


Inserting these into equation (Bl) gives 


p yj ( D w* + w >5). -P. S* ^ ±(e d 22 + (B2) 

\ 9x 8y 3z / ax 3y \ ay / 


Next the independent variables are formally transformed from x, y, and z to X, 
Y, and Z, where 
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Z = - 
L 


Thus, viewing U = U(X, Y, Z), one has 


dU=2HdX + £HdY + iI?dZ 
dX 3Y 0Z 


Hence, from equation (B3) it follows that 

au = 3U ax = i su 
ax ax ax l ax 

and similarly for the other derivatives. Then equation (B2) becomes 


p„ u i 


V ax ay az / l ax T 2 aY\ aY ) 


Finally, (B5) is made nondimensional by dividing every term by p^u^/L to obtain 


V 3X 3Y dZ 1 3X Re SY \ SY / 


+ . . . 


PorPl 


Pc o U « L 
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APPENDIX C 


RELATIONS INVOLVING MACH NUMBER M^ 

The hot gas inlet Mach number and the stagnation enthalpy at the inlet are by 
definition 


M 


oo 




(Cl) 


(C2) 


But for the case where the enthalpy per unit mass depends only on temperature 


v<.=/>p ( « < C3) 

1 00 

The power law representation of the relation between specific heat and temperature is, 
in terms of the dummy variable 


V« ) = C P,»(tJ 


l 


(C4) 


Performing the indicated operations and combining equations (C2) to (C4) give 


9 2 T_c_ m 

2 _ 00 p,°° 

” l+l 


. j 

V T J 


(C5) 


The local speed of sound c is 



Using the combined first and second laws of thermodynamics results in 

T dS = di - I dP 
P 


(C6) 


(C7) 
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But 



and differentiating the perfect gas equation of state gives 


(C8) 


dT = — - — dp (C9) 

R P p 2 R 

Inserting equations (C8) and (C9) into (C7) and setting dS = 0, one can solve for the right 
side of equation (C6) and find 


c 2 - C P» °° T °° 

'*'00 

C P,°° 1 

R 

Combining equations (Cl), (C6), and (CIO) yields 



Solving equation (Cl 1) for Tq/T^ gives, in implicit functional form, 


: = fi (“»’ - *) 


The Euler number is defined by 


Eu = 


Poo^i 


Using the perfect gas equation of state and equation (C5) in (C13) leads to 


(CIO) 


(Cll) 


(C12) 


(C13) 
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I + 1 


(C14) 


Eu =. 


2 

R 


' M . , 
_ 


Taking cognizance of equation (C12), one can write equation (C14), in implicit 
form, as 


E “ = * 2 (“«’ ^ l ) 


The Eckert number is defined by 


u 


Ec = 


c p,°° T o 


Combining equations (C5) and (C16) results in 


Ec =. 


l + lTn 


v l + l 


T J>) -1 

lvtJ . 


Hence, equation (C17), in implicit form, can be written as 


Ec = f 3 ( M co, ^ A 


(C15) 


(C16) 


(C17) 


(C18) 
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Figure 1. - Schematic drawing of air-cooled vane. 
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